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ABSTRACT 

We present new NuSTAR and Chandra observations of NGC 3393, a galaxy reported to host the 
smallest separation dual AGN resolved in the X-rays. While past results suggested a 150 pc separation 
dual AGN, three times deeper Chandra imaging, combined with adaptive optics and radio imaging 
suggest a single, heavily obscured, radio-bright AGN. Using VLA and VLBA data, we find an AGN 
with a two-sided jet rather than a dual AGN and that the hard X-ray, UV, optical, NIR, and radio 
emission are all from a single point source with a radius <0.2". We find that the previously reported 
dual AGN is most likely a spurious detection resulting from the low number of X-ray counts (<160) 
at 6-7 keV and Gaussian smoothing of the data on scales much smaller than the PSF (0.25" vs. 0.80" 
FWHM). We show that statistical noise in a single Chandra PSF generates spurious dual peaks of 
the same separation (0.55 ± 0.07" vs. 0.6") and flux ratio (39 ± 9% vs. 32% counts) as the purported 
dual AGN. With NuSTAR, we measure a Gompton-thick source (Vr = 2.2 ± 0.4 x 10^^ cm“^) with 
a large torus half-opening angle, 6>tor = which we postulate results from feedback from strong 

radio jets. This AGN shows a 2-10 keV intrinsic to observed flux ratio of ^150 (T 2 -iokeVint = 

2.6 ± 0.3 X lO^^ergs”^ vs. T 2 - 10 keVobserved = 1-7 ± 0.2 X lO^^ergs”^). Using simulations, we find 
that even the deepest Chandra observations would severely underestimate the intrinsic luminosity of 
NGC 3393 above z > 0.2, but would detect an unobscured AGN of this luminosity out to high redshift 
{z 5). 

Subject headings: galaxies: active — galaxies: Seyfert—X-rays:galaxies— galaxies:incliviclual (NGC 
3393) 
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1. INTRODUCTION 

The detection and measurement of the frequency of 
dual active galactic nuclei (AGN) is an important test 
of the merger-driven AG N model. Inspired by the work 
of Sanders et ah (1984), theorists invoke major galaxy 
mergers as the trigger for both substantial starbursts 
and quasar phases that sweep galaxies clear of gas, thus 
starvi ng both star formation a nd black holes of their fuel 


(e.g., Di Matteo et ah 2005). While secular processes 
rather than mergers dominate in lower-luminosity AGN 


(e.g., U~et al.|20Q6 |Schawinski et al.|2011 Treister et al 

20121), more powerful AGW phases have been seen tc 


2012), more powerful AG^ phases have been se en to 

occur in major, ga s-rich mergers in both near by ([Kq^ 


et al.| 2010 2011a) and high-redshift galaxies (Treister 


et al. II 2010). Currently, there is significant debate about 


whether mergers play a significant role in AGN activity 
particularly at high r edshift (e.g., Clsternas ^^ ^ ^ 


Kocevski et al. 
Alexander! 120151) . 


Schawinski et al.||2012[ [Brandt ^ 

Since AGN are rare, studies of dual 
AGW, closely separated galaxy nuclei (<30 kpc) both 
hosting AGN in an ongoing merger based on closeness in 
redshift and disturbed morphologies, have played an im¬ 
portant role in studying the frequency of merger-driven 
AGN activation. Over the last decade, dozens of these 
dual AGN on kpc scales have been found serendipitously 
in interacting galaxies using X-ray i maging, optical spec 
troscopy, or radio observations (e 


Barth et al. 1120081 jComerford et al. |2009 

|Koss et al. 

201 lb| Fu et al. 201 1|) or using lar 

ge-sca 

e surveys in 


However, despite ex tensive effort ( [iwasawa et al.^Oll 


Burke-Spolaor 2011) and simulations suggesting a peak 


of lu minous dual AGN activity at small separations (<10 
kpc. Van Wassenhove et al.|2012 ), few definitive sub-kpc 
dual AGIN have been found 
As the final phase in a major merger, many models 


suggest two heavily obs cured AGN should form (Hop 
kins & Hernquist||2009|). NGG 6240, the prototype of 
a close dual AGJN detected in the X-rays over 10 years 
ago, hosts two Gompton-t hick AGN separated by 1.5 


kpc (Komossa et al. 2003). The spectral shape of the 
10-100 keV band is especially useful for modeling such 
heavily obscured, Gompton-thick AGN because the 10- 
20 keV emission will decrease while higher energy emis¬ 
sion experiences less absorption when the AGN is in the 
Gompton-thick (>10^^ cm“^ ) to heavily Gompton-thick 


(5 X 10^^ cm“^ ) regime. With the new focusing optics 
on the Nuclear Spect roscopic Telescope Array {NuSTAR] 
[Harrison et al.|2013), the 10-79 keV energy range can be 
studied at sensitivities more than 100 x higher than pre¬ 
vious coded aperture mask telescopes such as the Swift 
Burst Alert Telescope (BAT) or INTEGRAL. 

NGG 3393 was reported to be a Gompton-thick dual 
AGN, with the smallest physical sepa ration (150 pc, 
0.6^0 ever observed using X-ray imaging (jFabbiano et al 
2011|). However, a recent study using adaptive optics 


(AU) on Subaru in the K-(2.2 /am) and V (3.8 /am) bands 
found no evidence of a secondary nucleus at 0.6" despite 
having higher reso lution (?^0.2") than Chandra (Iman- 
ishi & Saito 2014). Additionally, an analysis of tne nu¬ 


clear region of NGG 3393 with the VLT Imager and Spec¬ 
trometer for mid Infrared (VISIR) in four different N- 
band filters (8 — 13 /am) found only one compact source 


at a resolution of 0.3"(|Asmus et al. 2014). NGG 3393 
was observed by NuSTAR as part of a survey to observe 
the most obscured nearby AGN detected with the Swift 
BAT based on their spectra above 10 keV (Koss et al. in 
prep). 

The host galaxy of NGG 3393 is a face-on spiral galax; 


with a Seyfert 2 nucleus (Veron-Getty & Veron 


1986 


i 


The nucleu s shows polarized broad emission lines (|Kay 
et al.|200^ ). Radio images of the inner-kiloparcsec region 
reveal a cor e plus an apparent double-sided jet (jGooke 
et al. 2000) with a total extent of «700 pc. Tne ra¬ 
dio emission is surrounded by the S-shaped [O ill] emis¬ 
sion from the narrow line region (NLR) imaged with the 
Hubble Space Telesc ope Hubble Space Telescope {HST, 


Schmitt et al. 2001|). This is consistent with the jets 
creating denser regions of gas on their leading edges. A 
detailed examination of the NLR based on ground-based 
and space-based imaging and sp ectroscopy, as well as ra¬ 
dio observations, is presented by Gooke et al. (2000). The 
NGG 3393 nucleus is also a source of water maser emis¬ 
sion (Kondrat^^£F^I][20M 2008). EeppoSAA observa¬ 
tions in 1997 suggested columns of Nh = 3 x 10^^ cm“^ , 
but the large Fe Ka line equivalent width, high ratio of 
[O III] to soft X-ray flux, and the excess above 2 0 keV sug¬ 
gest a Gompton-thick AGN (jSalvati et al.|1997|). The col¬ 
umn was later found with XMM- Newton and BemoSAX 
data to be Nh = 4 x 10^"^ cm“^ (Guainazzi et al.||2005) 
implying that NGG 3393 is Gompton-thick, but not so 
obscured to suppress the nuclear emission fully. 

In this article, we test for the significance of a dual 
AGN as well as measure the strength of the AGN ob¬ 
scuration and reflection in NGG 3393. Since the original 
NGG 3393 study where a dual AGN was purportedly 
detected in the 6-7 keV range (Fabbiano et al. 2011), 
there are 340 ks of recently obtained Chandra observa- 
tions which increase the depth of 6-7 keV imaging by 
over a factor of three. In addition, there are recently 
obtained high-resolution milliarcsecond Very Long Base¬ 
line Array (VLBA) observations. Section 2 describes our 
imaging, spectra, and simulations as well as the analy¬ 
sis technique. We discuss the host galaxy optical and 
radio morphology in Section 3.1. Our results then fo¬ 
cus on using multi wavelength data to test whether the 
source is indeed an X-ray detected dual AGN (Section 

3.2) . We then discuss our results from fitting AGN mod¬ 
els to NGG 3393 and their implications for the geometry 
of the torus (Section 3.3). Finally, a summary and dis¬ 
cussion of the difficulty of detecting Gompton-thick AGN 
like NGG 3393 and its importance for higher redshift sur¬ 
veys is discussed (Section 4). Throughout this work, we 
adopt Ltrn= 0-27, ^a= 0.73, and Hq = 71 km s“^ Mpc“^. 
We assume a redshift for NGG 3393 of 0.0125 based on 
21-cm neutral hydrogen line measurements corresp ond- 
ing to a distance of 53 Mpc (Theureau et al.||1998 ). At 
this redshift, 1" subtends 250 pc. 

2. DATA AND REDUCTION 

Here we present an analysis of new imaging and grat¬ 
ing spectroscopy from the X-rays from NuSTAR (Section 
2.1) and Chandra (Section 2.2). We have also done a re¬ 
analysis of past X-ray observations of NGG 3393 (Section 

2.3) . A summary of the X-ray observations can be found 
in Table [2 Additionally, we have obtained near-infrared 
AO imaging (Section 2.4) and optical spectroscopy (Sec- 
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tion 2.5). Finally, we have also analyzed recently ob¬ 
tained VLBA observations and archival VLA radio ob¬ 
servations (Section 2.6). Errors are quoted at the 90% 
confidence level unless otherwise specified. 


Table 1 

X-ray Observations of NGC 3393 


^Effective exposure time (ks) after data cleaning and correction 
for vignetting. 

^Energy range studied in analysis. 


2.1. NuSTAR 

The raw data were reduced using the NuSTARDAS soft¬ 
ware package (version 1.3.1) jointly developed by the 
ASI Science Data Center (ASDC) and the California 
Institute of Technology. NuSTARDAS is distributed with 
the HEAsoft package by the NASA High Energy Astro¬ 
physics Archive Research Center (HEASARC). We ex¬ 
tracted the NuSTAR source and background spectra us¬ 
ing the nuproducts task with the appropriate response 
and ancillary files. Spectra were extracted from circu¬ 
lar regions 40" in radius, centered on the peak of the 
centroid of the point-source. The background spectra 
were extracted from three regions on the same detector 
as the source. The final exposure time after screening 
was 15.7 ks, with 821 and 758 background subtracted 
counts, respectively, for the EPMA and EPMB modules. 
The NuSTAR spectra are binned to a minimum of 20 
photons per bin using HEAsoft task grppha. 


Date 

Obsid 

Telescope 

Instrument 

Exp.^ 

(ks) 

Energy^ 

(keV) 

1997-01-08 

50035003 

BeppoSAX 

MEGS 

14.6 

0.5-10 

1997-01-08 

50035003 

BeppoSAX 

PDS 

14.6 

15-150 

2003-07-05 

140950601 

XMM 

PN, MOS 

15.8 

0.5-10 

2004-02-28 

4868 

Chandra 

ACIS-S 

29.7 

0.5-8 

2007-08-14 

702004010 

Suzaku 

PIN 

55.2 

15-40 

2011-03-12 

12290 

Chandra 

ACIS-S 

67.2 

0.5-8 

2012-02-29 

13967 

Chandra 

ACIS-S HEG 

179.8 

0.5-8 

2012-03-06 

14403 

Chandra 

ACIS-S HEG 

79.0 

0.5-8 ; 

2012-04-02 

14404 

Chandra 

ACIS-S HEG 

57.7 

0.5-8 ^ 

2012-04-08 

13968 

Chandra 

ACIS-S HEG 

28.5 

0.5-8 

2013-01-28 

61205001 

NuSTAR 

FPMA/B 

15.6 

3-70 ; 

2004-2010 

70 month 

Swift 

BAT 

7333 

14-195 p 



Figure 1. Radial profiles of the ACIS images (12290 and 4868) 
and zeroth order grating data (13967, 13968, 14403, 14404). A 
dashed line indicates the position of the Chandra hook asymmetry. 
We see no statistically significant evidence that any of the images 
are different from each other in terms of their PSF size. 


Eabbiano et al.||2011). The total number of counts at 6- 
7 keV in a 2" radius aperture increases from 145 to 449, 
a factor of 3.1, consistent with the increase in exposure 
time. We find no evidence that the imaging point spread 
function (PSE) is significantly worse in any of the ACIS 
images or zeroth order grating data at 6-7 keV (Eigure 

0 - 

The data were analyzed using Cl AO (Version 4.5) with 
the latest CALDB 4.5.8 provided by the Chandra X- 
ray Center (CXC). The subpixel positioning was applied 
using Ene rgy Dependent Snbpixel Event Reposi t ioning 
(EDSER, ^ 

with CIAO 


Tsunemi 
tool 


et al. 


2001 


Li et 


al 

TlTe 


2003 2004) 


acis_process_events. The events were 
screened for high-background periods with the CIAO tool 
deflare with a 2.5cr cutoff. Spectral responses were gen¬ 
erated using the CIAO tool specextract. 

In order to correct for offsets in the relative astrome¬ 
try of each Chandra observation, we use the CIAO tools 
fluximage to create exposure-corrected images and ex¬ 
posure maps, mkpsfmap to compute the point spread 
function (PSE) size across the image, and wavdetect 
for source detection. We include only the highest signifi¬ 
cance point sources when computing astrometry (SN>10, 
PSE size<1.5"), and detect seven sources in both images, 
excluding the nucleus. We find an offset of 0.41"±0.08" 


2.2. Chandra 

Archival Chandra Advanced CCD Imaging Spectrom¬ 
eter (ACIS) X-ray imaging taken on 2004 Eebruary 28 
(ObsID 4868 for 30 ks total exposure) and 2011 March 
12 (ObsID 12290 for 70 ks total exposure) is used in the 
analysis. The effective area of th e Chandra High Energy 


Transmission Grating (HETG, [C 


et al 


2000 ) 


zeroth-order image is about 60% of the ettective area of 
ACIS-S at 6.4 keV, so we also make use of 340 ks of recent 
HETG observations taken in 2012 (PI: Evans; ObsID 
13967, 170 ks; ObsID 13968, 30 ks; ObsID 14403, 79 ks; 
ObsID 14404, 61 ks), providing an equivalent of 210 ks 
ACIS-S imaging in the Ee Ka band («6-7 keV). There¬ 
fore, the combined depth of our imaging data reaches 
310 ks, which improves the exposure by a factor of 3.1 
in the 6-7 keV range where a dual AGN was purportedly 
detected in a past study of NGG 3393 (100 ks depth; 


in R.A. and 0.11"±0.05" in deck between ObsID 4868 
and 12290. We can also compute the relative astromet¬ 
ric offset using the 4-8 keV centroid of NGG 3393 and 
find similar values with the same offset direction (0.50" 
in R.A. and 0.07" in deck at 4-8 keV and 0.45" in R.A. 
and 0.12" in deck at 6-7 keV). This is consistent with the 
90% uncertainty (0.63") in ACIS-S astrometry of on-axis 
sourcej^ All seven sources show the same offset direc¬ 
tion in both R.A. and deck We then use the CIAO tool 
reproject_aspect to align images to the same world co¬ 
ordinate system (WCS). In images taken with the HETG 
offset sources are not available for crossmatching, so we 
combine the images based on the centroid of NGG 3393. 

NGC 3393 has evidence of extended emission, so we 
extract a nuclear region (2" radius) enclosing the PSE 

See http://cxc.harvard.edu/cal/ASPECT/celmon/ for more 
details on astrometry. 
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and a 40" radius matched to the NuSTAR observation. 
There is no evidence of pileup in any of the Chandra 
imaging or grating observations because of the low count 
rate of the nuclear emission 0.015 cts s“^in ACIS-S 
imaging). The nuclear region is corrected for the small 
amount of PSF flux (< 10%) beyond the extraction ra¬ 
dius. 

We also use the HETG Medium Energy Grating 
(MEG) and High Energy Grating (HEG) data from re¬ 
cently obtained archival data with a total exposure time 
of 340 ks. The gratings operate simultaneously, with the 
MEG/HEG dispersing a fraction of the incident photons 
from the two outer/inner High Resolution Mirror Assem¬ 
bly (HRMA) shells. The Ist-order HETG spectral prod¬ 
ucts were extracted using the tgcat pipeline for grating 
spectra in ISIS. We use a full-width of 4" in the cross¬ 
dispersion direction of NGG 3393 and extract it as a 
point source. This corresponds to a maximum distance 
from the central source of ?^500 pc. We find that the 
count rates in the MEG, HEG, and Ee K region for ob¬ 
servations 13967, 13968, and 14403 are all within la of 
each other. Observation 14404 has a larger difference, 
with a count ratio of l.Scr lower than the average of the 
others, which we do not consider significant. We there¬ 
fore coadd all the spectra. The grating spectra then have 
a total of 1259 counts in the MEG and 736 counts in the 
HEG and are fit with models simultaneously to improve 
the fit. We note that NGG 3393 is extended in the soft 
X-rays (< 3 keV) causing a degradation of the apparent 
spectral resolution. 

2.3. Other X-ray Telescopes 

We also include data from other X-ray observatories 
because of their overlap with NuSTAR and to provide 
constraints on variability between the energy range of 
0.5-200 keV. A summary of these additional X-ray ob¬ 
servations can be found in Table [TJ 

NGG 3393 was observed by BeppoSAX on 1997 Jan¬ 
uary 8 with the Low Energy Goncentrator Spectrometer 
(LEGS), the three Medium Energy Concentrator Spec¬ 
trometers (MEGS), and the Phoswich Detector System 
(PDS). There was no detection in the LEGS, so we use 
only the MEGS and PDS here. The MEGS contains three 
identical gas scintillation proportional counters, with an¬ 
gular resolution of ^0.7' full width half max (EWHM) 
and ^2.5' half power diameter (HPD). The MEGS event 
files were screened adopting standard pipeline selection 
parameters. Spectra were extracted from 4' radii aper¬ 
tures and the spectra from the three units were com¬ 
bined after renormalizing to the MECSl energy-PI rela¬ 
tion covering the energy range 0.5-10 keV. The PDS has 
no imaging capability to speak of, but has sensitivity to 
15-150 keV with some overlap with NuSTAR. The PDS 
data were calibrated and cleaned using the SAXDAS soft¬ 
ware, adopting the fixed Rise Time threshold method for 
background rejection. 

NGG 3393 was observed by Suzaku with the Hard X- 
ray Detector (HXD) PIN on 2007 August 14 with the 
nominal HXD pointing. We reprocessed the unfiltered 
events file using the standard Suzaku pipeline ETOOL 
aepipeline. The ETOOL hxdpinxbpi was used to ex¬ 
tract the source spectrum and to create a total back¬ 
ground spectrum from the tuned non-X-ray background 
for the NGG 3393 observation provided by the Suzaku 


team. The filtered events file, corrected for dead time, 
has a total exposure of 45 ks covering the energy range 
14-50 keV. The Suzaku source spectrum was binned such 
that, in each bin, the source is three times the back¬ 
ground. 

NGG 3393 was observed by XMM-Newton on 2003 July 
5. We processed XMM-Newton data using SAS (vl3.5.0) 
and the xmmextractor metatask for end to end process¬ 
ing of all the XMM-Newton raw data to science level 
data with filtering. This left 13.1 ks usable exposure. 
The background events are extracted from a annulus of 
inner radius 70" and outer radius 140" centered on the 
source position. We processed the EPIG MOS and PN 
data specifically using standard values for a point source. 
We use a 40" radius matched to the NuSTAR observa¬ 
tion for extraction and use the 0.5-10 keV range in EPIC 
MOS and PN data. There is no evidence of pileup in the 
observation. 

We finally use the Swift BAT stacked spectra and light 
curve from the publicly available 70 month catalog which 
spans 2004-2010 covering 15-195 keV. The light curve 
data was binned into 3 month intervals. Detail s of the 
BAT catalog reducti on procedure can be found in |Baum- 


gartner et al. (2013). 

When coniparmg observations between X-ray tele¬ 
scopes to assess variability it is important to consider the 
much wider PSEs of other instruments because of con¬ 
tamination by other sources as well as cross-calibration 
before comparing them with NuSTAR. This is es¬ 
pecially important for BeppoSAX data because of its 
lower angular resolution. Eor instance, the BeppoSAX 
MEGS data has a beam size of 4'. Based on XMM- 
Newton image, we find that other AGN separated be¬ 
tween 0.9' and 4' from NGG 3393 have ~45% of the 
4.5-12 keV counts. The BeppoSAX PDS EWHM is a 
large 1.3°. Based on the higher resolution Swift BAT 
maps, there is another source ?^0.88° NE detected at 
SNR=3.75 at R.A.=162.656° and decl=-24.306°, which 
is near (2.8' offset) the source IRXS J105027.3-241644. 
The source has a count rate roughly 29% of counts of 
NGG 3393 based on the Swift BAT 14-195 keV maps 
(0.000327 cts s“^ vs. 0.000796 cts s“^). 

We assume a cross-calibration of 0.93 between Suzaku 
PIN and NuSTAR and 1.25 between XMM-Newton and 
NuSTAR based on recent tests of other calibration 
sources. These values are based on the current best es¬ 
timates for NuSTAR (Madsen et ah, submitted). Eor 
BeppoSAX MEGS we assume a factor of 0.55 based on 
the contamination and a further cross-calibration of 0.8 
between MEGS and PDS following the instrument guide¬ 
lines. We have assumed a value of 1 for the cross¬ 
calibration of NuSTAR compared to Swift BAT because 
Swift BAT takes several years to reach the same signif¬ 
icance as NuSTARand has not yet been well calibrated 
(K. Madsen, 2015, Private Gommunication). 


2.4. High Resolution Optical and Near Infrared Imaging 

We imaged NGG 3393 on 2013 November 10 in the 
near infrared (NIR) using AO in laser guide star mode 
with the Near Infrared Camera 2 (NIRC2) instrument 
on the Keck-2 telescope. We used the wide camera with 
40 mas pixel”^ and a 40" field of view (EOV). We used 
a 3 point dither pattern for 18 minutes total in the Kp 
filter. Images were combined using a Strehl weighting 
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with Strehl ratios between 5-7% for a nearby star. The 
final PSF (FWHM) is 0.15", corresponding to physical 
scales of ^40 pc, as measured from nearby stellar sources 
in the final co-added image. The absolute astrometry of 
the NIR image was recalibrated using three bright stars 
from a PSl r-band image. 

We analyzed HST imaging from GO program 12185 
using Wide Field Camera 3 (WFC3) for filters F336W, 
F438W, and F814W. We also processed imaging from 
HST Advanced Camera for Surveys (ACS) F330W (Pro¬ 
gram 9379) and narrowband WFC3 FQ508N [O ill] (Pro¬ 
gram 12365). We apply AstroDrizzle to flat-field cal¬ 
ibrated images to create mask files for bad pixels and 
cosmic rays. Additionally, when the observations have 
multiple images, we drizzle-combine the input images us¬ 
ing the mask files, while applying geometric distortion 
corrections, to create a distortion-free combined image. 
The absolute astrometry of each image was finally recali¬ 
brated using three bright stars from a PSl r-band image. 


2.5. Optical Spectroscopy 

We also observed NGC 3393 using the University of 
Hawaii 2.2m telescope and the SuperNova Integral Field 
Spectrograph (SNIFS) on 2013 April 27 for a total du¬ 
ration of 3600 seconds. SNIFS is an optical integral field 
spectrograph with blue (3000-5200 A) and red (5200- 
9500 A) channels and a resolution of 360 kms“^. The 
SNIFS reduction pipeline SNURP was used for wavelength 
calibration, spectro-spatial fi at-fielding, cosmic ray re¬ 
moval, and flux calibration (Bacon et al. 2001 Alder- 
ing et al.||2006|). A sky image was taken after each 
source image and subtracted from each Integral Field 
Unit (IFU) observation. Flux corrections were applied 
based on the standard star Feige 34. We fit the spectra 
using an extensible spectroscopic analysis toolkit for as¬ 
tronomy, PYSPECKIT, which uses a Levenberg-Marquardt 
algorithm for fitting. 


2.6. Radio Observations 

NGC 3393 was observed with the VLA on 1992 Novem¬ 
ber 29 (project code AB618) at 1.4 GHz (3.0 x 1.5" beam 
size, pA= —14.1°), 4.9 GHz (0.8" x 0.4" beam size, 
PA= -7.7°), and 8.4 GHz (0.5" x 0.3", PA= -1.5°). 
In Figure we show 8.4 GHz contours overlaid on the 
color-scale 4.9 GHz VLA image (left panel). Three com¬ 
ponents, referred to as A (nucleus), B (NE jet lobe), and 
C (SE jet lobe) are clearly detected in agreem ent with 
past studies (See Section l,|Cooke et al.||2000). The po¬ 
sitions of the different components measured from the 
VLA 8.4 GHz map (Hie one with the highest resolution) 
are shown in Table The larger beam in the 1.4 GHz 
VLA map prevents us from resolving components A, B 
and C separately, and thus we no longer consider it in 
our analysis. 

We also reduced and combined two days of archival 
VLB A observations (2011 December 7-8; project code 
BS214) at 2.3 GHz (11.3 x 7.4mas, PA= 0.4°) and 
8.4 GHz (2.8 X 1.8mas, PA= 1.8°). Whilst A, B, and 
C are clearly detected in the VLA images (Eigure 1), 
only component B is detected in the VLB A image at 2.3 
GHz, and appears as an extended source (5'=2.3GHz = 
7.54 ±0.42 mJy). This further suggests B is a jet hotspot 
component. Deeper images are needed to detect compo¬ 


nents A and G. Alternatively, observations at higher fre¬ 
quencies can also help in the detection of component A, 
if this represents a case similar to 4G39.25, whose core 
has a spectr al turnover at high f requencies, between 43 
and 86 GHz ( Alberdi et al.|p^^ ). 


Table 2 

NGC 3393 radio components from VLA 8.4 GHz observations 


Gomponent 

Peak position 
RA(J2000) DEC(J2000) 

-^Jet —core 

(arcsec) 

A 

10 48 23.467 

-25 09 43.49 

- 

B 

10 48 23.406 

-25 09 44.06 

1.00 

G 

10 48 23.532 

-25 09 42.91 

1.06 


Position uncertainties in R.A. and decl. are better than ±10 mas. 


Table 3 

Flux densities and luminosities from VLA observations at 
matched resolutions. 


V 

(GHz) 

Gomponent 

Su 

(mJy) 

(102'^ erg s-i Hz-1) 

8.4 

A 

0.70 ± 0.04 

2.36 ± 0.14 


B 

9.22 ± 0.46 

30.98 ± 1.55 


C 

1.92 ± 0.10 

6.44 ± 0.33 

4.9 

A 

0.83 ± 0.05 

2.78 ± 0.17 


B 

15.16 ± 0.76 

50.95 ± 2.55 


G 

2.89 ± 0.15 

9.71 ± 0.49 


Uncertainties in flux density include a conservative 5% error in 
the point source calibration and the rms noise in the map, added 
in quadrature. 


In order to obtain useful spectral index information 
from the VLA images, we created 4.9 and 8.4 GHz maps 
with matched baseline ranges in wavelength, and the 
same convolving beam (0.65" x 0.35", PA= 0°). We 
then obtained a spectral index distribution map consid¬ 
ering only emission above 15a (see Table to avoid 
confusion between the different components (Eigure [^. 

3. RESULTS 

3.1. Single or Dual AGN? 

Before studying the X-ray spectra, we assess the signif¬ 
icance of the previously claimed dual AGN. Our results 
first focus on the high resolution radio, optical, and NIR 
morphology of NGG 3393 to test whether their is evi¬ 
dence of the purported dual AGN or a recent merger. We 
then test for astrometric offsets in the multi-wavelength 
data. Einally, we use deconvolution, modeling, and sim¬ 
ulations of Chandra X-ray data to test the significance 
of the purported dual AGN. 

3.1.1. Radio Morphology 

The structure recovered with the VLA is indicative 
of that from an AGN with a two-sided jet rather than 
a dual AGN. Gomponent A (nucleus) with a flat spec¬ 
trum is the core, whereas B and G have steep spectra 
with a flux difference easily explained by Doppler boost¬ 
ing. In this interpretation, components B and G are the 
approaching and receding components, respectively (see 
Eigure [^. Using the flux density enclosed by the 15cr 
contour m each component, we obtain spectral indices 
g=- 0.30±0.15, -0.90 ±0.13 and -0.74 ±0.13 for com¬ 
ponents A, B and G, respectively {Su oc z/^). 
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Figure 2. NGC 3393 VLA 8.4 GHz contours overlaid on color-scale 4.9 GHz map (left). Width of the radio image is 4.8". VLBA contour 
images from each radio component A, B and G at 2.3 and 8.4 GHz (right). The plus signs in the VLBA images represent the A, B and 
C peak positions as obtained from the VLA 8.4 GHz image, with an accuracy better than 10 mas. The circle in the VLBA images of 
component A, represents the dynamical location of the BH obtained from water maser emission ( [Kondratko et al.|20Q8| >. 



Figure 3. Pixel by pixel spectral index distribution between 4.9 
and 8.4 GHz emission detected with the VLA. The structure recov¬ 
ered with the VLA resembles an AGN with a two-sid ed jet rather 
than a dual AGN (both here and in|Gooke et al.|2000|). Gomponent 
A (nucleus), with a flat spectrum, is tlie core, and the steep spec¬ 
trum components are the lobes. Their flux difference explained by 
Doppler boosting, would be approaching and receding components, 
respectively. 


3.1.2. High-Resolution Optical and NIR Morphology 

We use the high-resolution HST imaging to search 
for signs of a recent merger. A HST tricolor (T336W, 
F438W, F814W) image can be found in Figure^ 

3393 is an early-type barred spiral t hat is nearly 


NGC 


lace on 


(SBars, [de Vaucouleurs et ^ 1995). NGC 3393 shows 


-25°09'00.0"^^^| 
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Figure 4. Tricolor HST images (F336W, F438W, F814W) of 
NGG 3393. The image is 2' on a side. NGG 3393 is a spiral 
galaxy that shows faint ti ghtly wound spiral arms which are much 
stronger in the Ha image (Gooke et al.|2000|) that connect to a bar 
at PA Psl60°. There is no evidence ot any asymmetry in the spiral 
arms as would be expected from a recent merger. 


faint, tightly wound, blue spiral arms, which are much 
stronger in the Ha image (Cooke et al. 2000). There is 
no evidence of any asymmetry in the spiral arms or any 
sign of tidal tails indicative of a recent merger. 

We then compare the HST images of the nuclear region 
of NGC 3393 to radio contours for the F330W, [O Hi], and 
F814W filters (Figure!^. In the nucleus, an S-shaped 
source can be seen in the bluer images coincident with 
the radio emission. In the UV and [Olll] there is ev¬ 
idence of a small biconical region on the nucleus that 
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is extended by 0.2". The identical biconical structure 
in the NUV and [O ill] suggests that the near-UV emis¬ 
sion is produced in the same region as the ionized gas, 
and thus the nebular continuum and the [NeV] emission 
must be the main contribu tors as noted by past studies 
(Munoz Marin et al. 2009). Additionally, the biconical 
structure in the HST W and [O ill] is likely due to star 
formation or photoionized emission because of their lack 
of detection in the NIR imaging. Finally, in the NIRC2 
images (Figure [^, we find a bar at PA?^160°, with a 
smooth distribution and no evidence of a secondary nu¬ 
cleus. 


3.1.3. Astrometry 

In a past study (Fabbiano et al.|2011), the optical HST 
image and water maser seen in the radio were found 
to be offset from each other, suggesting the presence 
of two AGN. We therefore compare the centroids of all 
the radio, optical, NIR, and X-ray imaging (Figure]^. 
The absolute position of the dynamical center was esti¬ 
mated to within 1 mas based on the mean position of 
the low-velocity maser f eatures (R.A.bh = 162 .09777°, 


decl.BH = —25.162077°; Kondratko et al. 2008) and we 


use this observation as the absolute reference. We find 
the hard X-ray, UV, optical, NIR, and UV emission are 
all coming from a single point source within a radius of 
< 0.2" (Icr). This is fully consistent with astrometric 
errors in the X-rays of ±0.3" (la) based on Sgr A* field 
and ±0.15" (la) in the Hubble Deep Field as well as well 
calibrated sky surveys such as 2MASS (0.3" accuracy for 
extended sources like NGC 3393 at larH 


3.1.4. Testing Chandra Imaging for a Dual AGN 

We next test for the significance of the two reported nu¬ 
clei using modeling and other techniques. The 0.5" pixel 
size of Chandra AGIS provides a limit to the resolution at 
which structures can be unambiguously detected. Super¬ 
resolution can be achieved using the high cadence aspect 
solution and binning sky pixels at sub-integer sizes which 
will, on average, give a good approximation to the actual 
arrival direction of the photon. However, the small sep¬ 
aration of the two proposed nuclei in past studies (0.6") 
is very close to the telescope resolution (0.6", 90% confi¬ 
dence limits using AGIS-S). Also, there are known PSF 
asymmetries that appear at small separations (<0.8"). 
Additionally, the previously claimed detection of a dual 
AGN was based on the 6-7 keV emission, where the col¬ 
lecting area of Chandra is low and the PSF is somewhat 
worse. Finally, the previously reported detection com¬ 
bined two images and a possible incorrect offset between 
the two images would appear as two distinct sonrcej^ 
Due to these complicating issues it is important to study 
the Chandra imaging in detail. 

We first study the native Chandra PSF with EDSER. 
To generate a simulate d PSE, we use the Chandra ray¬ 
tracing program ChaRT (jGarter et al.|2003|) and the MARX 
software version 5 to project the ray-tracings onto the 
AGIS-S detector. ChaRT takes as inputs the position of 


For details on the accuracy of astrometry with Chandra, see 
http: //cxc .harvard. edu/cal/ASPECT/improve_astrometry.html. 

For more details on the Chandra PSF and its asymmetry, see 
http: //cxc .harvard. edu/cal/Hrc/PSF/acis_psf _2010oct .html 
and http: //cxc. harvard. edu/ciao/caveats/psf .artifact. html. 


the point source on the chip, the exposure time, and 
the spectrum of NGG 3393. To create simulated images 
of a spatially offset dual AGN, we use the GIAO tool 
reproject_events. The Chandra PSE was created by 
simulating a point source at the same position on the 
detector in each observation separately. We run an on- 
axis ray trace simulation at 6.4 keV with 20000 counts 
using pixadj=EDSER and AspectBlur=0.19", which rep¬ 
resents the telescope pointing uncertainty. We find the 
Chandra PSE to have a EWHM of 0.8"when fit by a Mof¬ 
fat, Gaussian, or radial profile. The radius of the 50% 
energy enclosed function is 0.5". 

To detect a dual AGN we look at the stacked as well 
as the individual images deconvolved by the expected 
PSE. The stacked image from 3-6 keV and 6-7 keV can 
be found in Eigurej^ With more than a factor of three 
deeper Chandra imaging at 6-7 keV the summed imaged 
shows no evidence of a secondary AGN either at 6-7 keV 
or 3-6 keV. The brightest source falls between the two 
previous detections and is consistent with a single AGN. 
We also investigate all of the images separately at 6-7 
keV and find no evidence of a dual AGN (Eigurej^. 

The deconvolved image for the deepest 70 ks exposure 
from the previously claimed dual AGN paper in differ¬ 
ent energy ranges can be found in Eigure We find 
extended emission in the 0.5-3 keV data (source extent 
1.43" ± 0.11" with PSE source extent of 0.47" ± 0.05"), 
but no significant extended emission is found from 3- 
6 keV (1.01" ± 0.18" for observation and 0.86" ± 0.13" 
for the PSE) and 6-7 keV (0.94" ±0.25" for nucleus and 
0.90" ±0.25" for PSE). 

We then studied whether the purported dual AGN 
could be generated by Gaussian statistical noise in a sin¬ 
gle PSE (Eigure with the same exposure as the 70 
ks detection. We ran nine simulations of a single simu¬ 
lated Chandra PSE. We used a Gaussian 2 pixel smooth 
kernel at 1/4 pixel binning following the sub-pixel bin¬ 
ning and smoot hing scales of the purp orted dual AGN 
discovery paper (jEabbiano et al.||2011). We inspect the 
smoothed Chandra image and locate the two brightest 


peaks, when there is an increase in the smoothed counts 
compared to adjacent pixels. We only include bright 
secondary peaks where the sum of the counts in a 
0.3" radius circle is greater than 25% of the sum 
of the brightest peak. Most of the images (6/9 or 
67%) show two bright peaks at random orientations de¬ 
spite being statistical noise from a single modeled PSE. 
This is because of the low counts (1-3 per pixel) and the 
fact that the smoothing kernel (0.25") is significantly 
smaller than the resolution at 6.4 keV (0.8"), causing 
Possion noise to appear as two distinct peaks. We find 
that the average separation of the two spurious bright 
peaks is 0.55" ± 0.07", with the secondary having on av¬ 
erage 39 ± 9% of the total counts. This is similar to the 
claimed dual detection in NGG 3393 in separation (0.6") 
and in the secondary counts (32% of total), suggesting 
the detection is statistical noise. 


3.2. Broadband Observations of a Compton-thick AGN 
3.2.1. High-Resolution Chandra Imaging and Speetroseopy 

We first examine the high resolution imaging from 
Chandra to determine the nature of soft extended emis¬ 
sion in order to best model our NuSTAR data. A plot 
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RA (J2000) 




Figure 5. HST image of the nuclear region of NGC 3393 compared to VLA 8.4 GHz radio contours for the 330 W, [Olll1. and 814W filt ers. 
A green circle of 0.2" radius represents the dynamical location of the BH obtained from water maser emission (|Kondratko et al.|2008|. In 
the nucleus an S-shaped source can be seen in the bl uer images coincide nt with the radio emission. The radio lobes created denser regions 
of gas on their leading edges similar to a bow-shock (|Gooke et al.|2000|). In the F330W UV and [Olll] there is evidence of small biconical 
region of 0.2" that is not detected in the redder bands. i'Jiere are sorne small absolute offsets with the radio data, however these offsets 
are consistent within the 0.2" astrometric uncertainty. 


of the nuclear PSF emission compared to a 40" region 
can be found in Figure Below 3 keV the majority 
of emission is extended (M%) based on comparing the 
background-subtracted count rates from a 40" region 
(0.056 ± 0.003 cts s“^) compared to a region enclosing 
the PSF (0.011 ± 0.001 cts s“^). In the 3-6 keV range, a 
significant amount of emission is still extended beyond 
the PSF emission (44%; 0.0034 ± 0.0008 cts s“^ com¬ 
pared to 0.0019 ± 0.0003 cts s“^). Between 6-8 keV, 
there is no evidence that any of the emission is 
extended (0.0019 ± 0.0008 cts s“^ compared to 
0.0014 ± 0.0003 cts s“^), consistent with our 2D 
MARX simulations. A steep soft power law component 
with Ts = 3.3 is a good fit to the extended Chandra data. 

The extended emission shows a significant peak at 
1.8 keV which suggests strong photoionized emission. 
We next focus on studying this emission at higher reso¬ 
lution using the grating spectroscopy from Chandra. We 
note that the grating spectra are from a region of only 
4" in the cross-dispersion direction, so they correspond 


to a region similar in size to the nuclear region discussed 
in the last paragraph, rather than from the 40" region 
where much of the extended emission below 3 keV 
is observed. Using the other ACIS-S imaging data, 
we find that this 4" region includes 53% of the total 
background subtracted counts found in a 40" region at 
0.5-8 keV. Gaussian line fitting was carried out using 
an automated pro cedure which fits the spectrum in 2A 
intervals following Kallman et al. (2014). The Gaussians 
are added to the model at wavelengths corresponding to 
known lines. A summary of detected emission lines can 
be found in Table IH 


The most prominent observed spectral features are an 
emission line at 6.39 keV with an equivalent width of 
1.11 keV associated with neutral iron, an emission line 
at 6.70 keV indicative of He-like iron with an equivalent 
width of 0.70 keV, and an emission line near 7 keV with 
an equivalent width of 1.13 keV which may be associated 
with a combination of the K absorption edge from near¬ 
neutral gas, emission from H-like iron (6.97 keV), and Fe 
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Figure 6. Top: Image of NGC 3393 in the K' band. The PSF 
FWHM is 0.10'', corresponding to physical scales of 25 pc. The 
images show a bar at PA 160°. Bottom: Zoom in around NGC 
3393. The diameter of the smallest white contour is ^^r^lOO pc. We 
do not see any evidence of an obscured secondary nucleus and the 
contours show a smooth potential consistent with a single nucleus. 


K/5 (7.06 keV). The high value of the equivalent width 
of the iron lines is consistent with Compton-thick AGN 


(e.g. [Krolik fc Kallman|p^^ Levenson et al.||2QQ2 ). 


3.2.2. NuSTAR Spectral Fitting 

We then examine the NuSTAR data (3-70 keV) alone. 
The NuSTAR FPMA and FPMB cross-normalizations 
were allowed to vary independently in the fit and are 
consistent with the 5% error from measured calibration 
sources. A simple absorbed power-law model does not 
fit the data well (x^/z^ = 121/72), but reveals that 
the spectrum is highly absorbed (A^h ^ 10^"^ cm“^ ) 
and hard (effective photon index ^ 1.2) - indicating 
Compton-thick absorption and a significant contribu¬ 
tion from a scattered (reflected) component. As a next 
step we use the appro ximate phenomenological models 
pi cabs (Yaq oob 1997|) for the transmitte d component, 
and pexrav (Magdziarz & Zdziarski||1995 ) to model the 
reflected component. T'he latter is implemented so that 
it reproduces only the reflection, R < 0 in Xspec termi- 
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Figure 7. Offsets in position of the galaxy nucleus in the X- 
ray, optical, and NIR as relative to radio position. We find no 
significant offset in any of the observations, but very small offsets 
in absolute astrometry (<0.2"), suggesting the emission is coming 
from a single AGN down to these scales (50 pc). The largest scatter 
is from the X-ray centroids, but a combination of all the X-ray 
observations to date finds a smaller scatter from the radio position. 
Additionally, when combining optical with X-ray imaging of offset 
sources, the X-ray observations show reduced offsets (e.g. 12290 
<0.1"). The optical (PSl) and 2MASS images also have offsets 
in absolute astrometry, but these can be reduced to < 0.05" by 
combining imaging of offset stars. 


Table 4 

NGG 3393 MEG and HETG Lines 


Line 

Observed 

(keV) 

Lab 

(keV) 

a 

(keV) 

Norm 

(lE-07) 

O VIII 

0.653T0.002 

0.654 

0.0010 

102.0 

Fe XVII 

0.727T0.004 

0.720 

0.0003 

16.5 

Fe XVII 

0.8I3T0.004 

0.826 

0.0007 

13.9 

Fe XVIII 

0.873T0.003 

0.873 

0.0009 

10.0 

Ni XIX 

0.884T0.003 

0.883 

0.0003 

5.9 

Ne IX 

0.906T0.003 

0.905 

0.0005 

6.5 

Fe XIX 

0.9I6T0.003 

0.917 

0.0004 

16.6 

Fe XIX 

0.924T0.003 

0.922 

0.0007 

13.6 

Ne X 

I.022T0.003 

1.022 

0.0015 

16.8 

Fe XXIII 

I.I3IT0.003 

I.I29 

0.0001 

3.5 

Fe XXIV 

I.I77T0.003 

I.I68 

0.0017 

3.6 

Mg XI 

I.330T0.003 

I.33I 

0.0003 

8.5 

Mg XI 

I.352T0.003 

1.352 

0.0039 

11.9 

Mg XII 

I.478T0.005 

1.472 

0.0074 

7.3 

Mg XII 

I.74IT0.003 

1.745 

0.0103 

8.6 

Si XIII 

I.837T0.003 

1.839 

0.0061 

2.2 

Si XIII 

I.86IT0.003 

1.865 

0.0078 

5.5 

Si XIV 

2.006T0.004 

2.005 

0.0060 

4.9 

Fe Ka 

6.388T0.03 

6.40 

0.0541 

43.4 

Fe XXV Ka 

6.70IT0.04 

6.70 

0.0519 

25.7 

Fe XXVI Ka 

7.004T0.05 

6.96 

0.0727 

38.1 


nology. We link the photon index of the two components 
and flx the inclination at cos i = 0.45 and the high-energy 
cut -off at 200 keV (following recent cut-off measurement s 


e.g. Gilli et al. 2007 Malizia et al. 2014 Ballantyne 2014). 
In order to account for the clear line-like excess around 
the energy of the neutral iron Ka line at 6.4 as well as 
the other iron emission lines at 6.7 and 7.0 keV we add a 
Gaussian component (cr = 10“^ keV) fixed at those en¬ 
ergies. The best fit using this model (x^/^ = 51/70) is 
found for F = 1.9T0.3 and Nn = (3.1T0.6) x 10^^ cm“^ . 

The model components used in the previous paragraph 
are only approximate; with the quality of the NuSTAR 
hard X-ray data we are able to also test more physi- 
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Figure 8. Stacked X-ray image with 3.lx deeper imaging rebinned to 1/4 pixel (0.125") resolution with 0.25" smoothing (top) or no 
smoothing (bottom). Highlighted as circles are the 6-7 keV centroid (purple), PSF FWHM (green), as well as the positions of the reported 
dual AGN from in past studies (blue and red). The red circle indicates the secondary AGN reported at 6-7 keV but not at 3-6 keV in past 
observations. The 3-6 keV image is on a higher linear scale because of the larger number of counts. We see no evidence of a dual AGN 
with the improved astrometry and deeper imaging, and the emission is consistent with the same single point source at both 6-7 keV and 
3-6 keV. 


cally motivated models of the Comp ton-thick obscuring 
torus . One such model is MYtorus (Murphy & Yaqoob 
2009). It includes both the transmitted and the scat¬ 
tered components, as well as fluorescent line emission 
computed for a torus with a half-opening angle of 60° 
through Monte Carlo simulations. A steep soft power 
law component with Tg = 3.3 is included in the model in 
order to account for a part of the complex soft emission, 
as MYtorus predicts negligible flux below 5 keV for any 
line of sight through a Compton-thick torus. The best fit 
for the NuSTAR data is found for the MYtorus model in 
the coupled modj^ with nearly edge-on inclination. A 

This is the default mode, in which the components of the 


formally good fit < 1) is found for fixed T = 1.9 

and any inclination greater than 70°. Fixing the inclina¬ 
tion at i = 80° we find = (2.5 ± 0.2) x 10^^ cm“^ for 
fixed r = 1.9, and = (2.1 ± 0.5) x 10^^ cm“^ for the 
fitted r = 1.7 ±0.2. 

Another physically motivated model for a torus, al¬ 
though the geometr y is not exactly toroidal in this 
case, is the BNtorus (jBrightman & Nandra|2011|) shown 
in Figure We again hx the inclination to nearly 
edge-on, at i = 87°. We find a statistically good fit 

model are linked so that they represent reprocessing of radiation by 
a uniform toroidal distribution of gas. For details on the different 
usage modes of the MYtorus model and their interpretation we refer 
the reader to Yagoobj ( |2012| ). 
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Figure 9. From left to right: data, gaussian smoothed, csmooth using Poisson statistic, and deconvolved image in the 6-7 keV range at 
1/4 sub pixel binning (0.125'') for different observations of NGC 3393. A green and red circle indicate the 50% and 90% energy encircled 
function. Note how the 3 pixel Gaussian smoothing shows features that are not present in the csmooth which accounts for Poisson statistics 
in the low count regime. The green shape indicates the position of a hook that is a known HRMA PSF asymmetry found in brighter 
observations that depends on the roll angle of the spacecraft plotted with GIAO tool make.psf .asymmetry.region. We find that there is 
excess emission offset from the PSF by 0.8" at 6-7 keV, but its orientation, separation, and brightness level (<10% of central PSF or 1-8 
counts) are consistent with the hook asymmetry feature of the HRMA found in brighter images. 
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Obs 12290 

lEnergy: 0.5 - 3 keV, Cts: 2533 
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Figure 10. From left to right: data, projected PSF, and deconvolved image for obsid 12290 using the predicted PSF of NGC 3393. The 
first row is 0.5-3 keV, the second row is 3-6 keV, and the bottom row is 6-7 keV. We find extended emission in the 0.5-3 keV data, but no 
significant extended emission is found at 3-6 keV or at 6-7 keV. The green shape indicates the position of a hook that is a known HRMA 
PSF asymmetry that depends on the roll angle of the spacecraft. The extended emission below 3 keV is coincident with the [O ill] and UV 
emission coming from the denser gas that is swept up by the leading edge of the radio jet. The low counts (110 counts) limit the ability to 
find a dual AGN at scales below '^1". 
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Figure 11. Nine simulations of a single simulated Chandra PSF shown unsmoothed and smoothed. We have used a Gaussian 2 pixel 
smooth kernel at 1/4 pixel binning following the sub-pixel binning and smoothing scales of in the purported dual AGN discovery paper 
(|Fabbiano et al.|2011|). Most of the images (6/9 or 67%) show two bright peaks, despite being statistical noise from a single PSF. The green 
circles are centered on the two brightest peaks in each image. This is because of the low counts (1-3 per pixel) and that the smoothing 
kernel (0.25") is significantly smaller than the resolution at 6-7 keV 0.8") causing Possion scatter noise to generate two distinct peaks. 
We find that the average separation of the two fake peaks is 0.55" ±0.07" with the secondary having on average 39± 10% of the total counts 
which is similar to the purported dual detection in NGG 3393 in separation (0.6") and the secondary counts (32% of total), suggesting the 
detection is statistical noise. 



Figure 12. Gount rate from a 40" aperture in Chandra (red filled 
squares), the same as the aperture used for NuSTAR and XMM- 
Newton^ compared to the central 2" (black circles). Below 3 keV 
the majority of the emission is extended, with a significant peak 
at 1.8 keV suggesting significant photoionized emission. In the 
3-6 keV energy range, a smaller amount of emission is extended 
beyond the PSF emission (44%). In the 6-8 keV energy range, there 
is no evidence that any of the emission is extended. Because of this 
contamination, we focus on the > 3 keV emission for measurement 
of AGN emission models. 


{x^l^ = 50/70) for the NuSTAR data with F = 1.8±0.2, 


A^h = (2.2 ± 0.4) X 10^^ cm ^ , and a half-opening 
angle of the torus 6>tor = degrees. The small 

upper limit is present because the spectrum changes 


rapidly at very wide opening angles (Brightman & Nan- 


dra 2011). The model yields an unabsorbed luminosity 
= 2.6 ± 0.3 X 10^^ 


erg s ^ 

observed I/ 2 -iokeV = 1-7 ± 0.2 x 10^^ ergs 


of ± 2 - 


lOkeV 


compared to an 
. The ratio 


of the unabsorbed to absorbed flux is 149. 

In summary, NuSTAR data are consistent with the pic¬ 
ture conveyed by both the MYtorus and the BNtorus 


models: a torus with relatively small covering factor 
(wide half-opening angle) seen through its Compton- 
thick wall with A^h ~ 2 x 10^^ cm“^ . Both the transmit¬ 
ted and the scattered components are needed for a good 
fit with nearly equal contributions above ^2 keV. We 
therefore conclude that the hard X-ray data clearly point 
towards an edge-on, borderline Compton-thick torus. 

The origin and the correct physical model for the soft 
X-ray emission are difficult to disentangle, since con¬ 
tributions from photoionization by the AGN and star 
formation in the galaxy, as well as Thomson-scattered 
AGN continuum, contribute a different fraction on dif¬ 
ferent spatial scales. Under the assumption that star 
formation alone is responsible for the infrared flux, the 
IRAS luminosity yields a l ow star formation rate of ~4 
M^/yr on 30 kpc scales (Kondratko et al. 2008), cor¬ 
responding to approximately I/ 2 -iokeV = lo^^ ergs“^, 
given the rela tionship between sty formation and X- 
ray emission (Lehmer et al. 2010). As this value is 
more than an order of magnitude below the observed 
A2-iokeV = 1-7 ± 0.2 X lO^^ergs”^, this suggests star 
formation is not an important contributor above 2 keV. 
The exact contribution of scattering, photoionization, 
and transmitted continuum to the soft X-rays is difficult 
to constrain because of model degeneracies other than 
that the photoionization lines dominate the emission be¬ 
low 2 keV based on the Chandra grating data. However, 
we find that different parameterizations for the soft part 
of the spectrum fit jointly between NuSTAR and Chan¬ 
dra lead to self-consistent measurements of the intrinsic 
photon index and the column density, which match the 
values found based on the NuSTAR data alone within 
their 90% uncertainties. 

In summary, we firmly establish that the hard X-ray 
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Figure 13. Fits to 340 ks of Chandra grating data from the cen¬ 
tral 2" aperture. Even in the 2" nuclear region, the majority of the 
emission between 1 and 2 keV is from photionized line emission. 
Top: MEG data (black data, red model) between 0.5-1.2 keV (10- 
22 A). Middle: Combined fit from Chandra grating spectra from 
MEG (black data, red model) and HEG (green data, blue model) 
between 1.2-2.5 keV (5-10 A). Bottom: Fit from Chandra grating 
spectra from HEG between 4.0-8.0 keV (1.5-3 A). 


data support a Compton-thick torus seen edge-on. This 
geometry agrees well with the fact that a disk water 
mega-maser has been detected in NGC 3393, since its 
detec tion implies a nearly e dge-on view of the masing 
disk ( Kondratko et al.||2QQ8 ). As the torus and the mas¬ 
ing disk may correspond to different physical scales, the 
hard X-ray modeling implies a high degree of alignment 
of the AGN sub-structures. 


3.2.3. Variability and Spectral Fitting with Additional X-ray 

Observatories 

We consider whether NGG 3393 is variable, consistent 
with a transmission-dominated AGN, or constant, as ex¬ 
pected for a reflection-dominated AGN. We test variabil¬ 
ity first using similar telescopes to avoid cross-calibration 
issues. We explored longer term variability in the hard 
X-ray band using the Swift BAT 70-month data taken be¬ 
tween 2004 and 2010. We find the full 14-100 keV band 
light curve of a region centered on NGG 3393, binned in 
3 month intervals, is consistent with a constant source 
model based on a chi-squared test (Figure f^. We can 
also exclude variability at levels greater a factor of two 
when the data is binned into year long timescales at the 
90% level. Splitting the Swift BAT data into several 
energy bands shows no significant variability (> 90%) 
either. In the Chandra data taken in 2004 and 2011- 
2012 we also find no variability within the observation 
or between the observations at the 90% level. This is 



5 10 20 50 

Energy (keV) 


Figure 14. Fit to NuSTAR data for NGG 3393 with BNtorus 
model. We fix the inclination to nearly edge-on, at z = 87° and find 
F = 1.8 ±0.2, A^h = (2.2 ±0.4) x 10^^ cm“^ , and the half-opening 
angle of the torus, ^tor = degrees. Similar properties within 

error are found from the MYtorus model as well. 

consistent with the findings of the Chandra Source Gat- 
alog which finds there is definitely no variability in the 
soft and medium energy bands (<2 keV) and likely no 
variability in the hard band based on a Gregory Loredo 
variability algorithm. Finally, within the XMM observa¬ 
tion we see no evidence of variability in the PN or MOS 
cameras at the 90% level which agrees with the 3XMM- 
DR4 XMM catalog. In summary there is no evidence 
of variability within any of the specific observations for 
NGG 3393. 

We finally fit all of the X-ray data using the BNtorus 
model and allow the normalization to float (Figure [T^. 
The data seem consistent with a constant flux between 
2003-2012, followed by a 3cr increase in flux when the 
NuSTAR observation was taken in 2013. INTEGRAL 
and Suzaku HXD have found Swift BAT cross-calibration 
to be systematically lower than their nominal values by 
the same factor of 0.82±0.03 factor at 14-100 and 14- 


70 ke V, respectively, using different studies (Molina et al. 
2013|). Given this systematic offset, the variability is 


only significant at a level of 1.8cr. There seems to be 
a bright state during the 1997 BeppoSAX PDS observa¬ 
tion, though it is difficult to exclude source contamina¬ 
tion in this observation because of the nearby source in 
the BAT maps. Further studies with the high sensitiv¬ 
ity of NuSTAR would be necessary to confirm this high 
energy variability and to better understand the cross¬ 
calibration using very bright sources. 

We finally refit all of the X-ray data using the BNtorus 
model to confirm our conclusions about the torus based 
on the NuSTAR data alone. We exclude the BeppoSAX 
data because of possible source contamination (see 2.3). 
With an inclination to nearly edge-on, at i = 87°, we 
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find a statistically good fit jv = 145/138) with F = 
1.82±0.09, A^h = (2.34±0.19) x 10^^ cm-^ , and the half¬ 
opening angle of the torus, ^tor = degrees. These 

values are all consistent within the errors to the NuSTAR 
only fits. The data is consistent with no high energy 
cutoff with a lower limit of Ec >169 keV when fitting a 
Pexrav and a high energy cutoff. 


3.2.4. Intrinsic Luminosity, Eddington Ratio, and X-ray 
Radio Loudness 

We can calculate the expected bolometric luminosity 
(Lboi) of NGC 3393 using optical and MIR measure¬ 
ments that can be compared to the X-rays. This is crit¬ 
ical because of the difficulty estimating the intrinsic lu¬ 
minosity based on the X-rays in Compton-thick AGN. 
We calculate an extinction-corrected [Olll] luminosity 
of of 1.5 X 10^^ ergs ^ from a 2.4" diameter region or 
3.9 X 10^1 erg s~ ^ from a 4.8" diamet er aperture from 
the SNIFS IFU. iLamastra et al.| (|2009|) reported the lu¬ 
minosity dependent bolometric correction factors. At 
^[oiii] = ~ lO^^ergs”^, the correction is 142, sug¬ 

gesting a bolometric luminosity of 2.1 x lO^^ergs”^ for 
a 2.4" diameter region or 5.5 x lO^^ergs”^ for a 4.8" 
diameter region which includes the extended [O ill] emis¬ 
sion. We use the extracted low resolution Spitzer spectra 
from the Cornell Atlas o f Spitzer IRS Sources (CASSIS; 
Lebouteiller et al.||2011), with a tapered column extrac- 
tion. We measure [JNev] and [Oiv] lines with luminosi¬ 
ties of 9.2 X lO^^ergs”^ and 5.7 x lO^^ergs”^, respec¬ 
tively. The [Ne v] luminosity corresponds to a bolometric 
luminosity of 8.3 x 10^^ erg s“^. Based on the unabsorbed 
2-10 keV luminosity of 2.6 x 10^^ erg s~^, we use the bolo - 
metric correction of 30 from |Vasudevan fc Fabiaii ( 2009[) . 
We find a bolometric luminosity of 7.8 x 10^^ ergs“\ 
which is consistent with the [Nev] estimates. 

We modeled the spec tral energy d i stribu tion (SED) 
of NGC 3393 using the Assef et al. (2010) 0.03-30/im 
empirical AGN and galaxy templates to understand the 
strength of the AGN emission in the IR. Each SED is 
modeled as a best-fit, non-negative combination of an 
old stellar component, a starburst, plus an AGN compo¬ 
nent. Only the AG N component is fit for d ust redden¬ 
ing (for details, see Assef et al. 2008, 2010). The red¬ 
dening corrected 6 fim luminosity of the best-fit AGN 
component is (1.5 ± 0.2) x lO^^ergs”^. The model¬ 
ing outputs I/6ytim, the derived intrinsic luminosity of 
the AGN component at rest-frame 6 jam, as well as the 
reddening of the AGN component, E{B — E)agn- Eor 
the ty pical gas-to-dust ratio observed by |Maiolino et ah 
(2001) for luminous AGN, the nuclear reddening values 
of e(B — V)agn = 12.6 ± 6.2 implies a gas column of 
Nh-1 X 10^^ cm^, consistent with the Compton-thick 
AGN interpretation. 

Combined with the measured mass of the SMBH har¬ 
bored by the AGN we can estimate the Eddington frac¬ 
tion, L^oi !where L^dd is the Eddington luminosity. 
VLBI maps of nuclear water maser emission have been 
used to measure the BH mass to be (3.1 ± 0.2) x 10^ Mq 
(K ondratko et al. 2008) with a corresponding Edding¬ 
ton luminosity of 4 x lO^^ergs”^. The Eddington ratio 
is then ^0.2 based on the [Nev] or intrinsic 2-10 keV 
emission, or ^0.1 based on the [Om] data. 

Based on the intrinsic 2-10 keV X-ray luminosity and 


the nuclear 4.9 GHz emission from component A, we 
can estimate the X-ray radio loudness parameter, Rx = 
iyLjy{b GHz)/Lx- Radio quiet objects are typ ically found 
at logRx < —4.5 ( iTerashinia fc Wilson|2o53 ). We find a 
value of -6.28 for NGG 3393 consistent with a radio quiet 
rather than radio loud object. 


4. SUMMARY AND DISCUSSION 

We study the Compton-thick AGN NGC 3393 using 
recently obtained NuSTAR, Chandra, VLB A, and Keck 
data combined with archival HST and VLA data. This 
nucleus has been claimed to harbor a dual AGN based 
on previous Chandra imaging. We find: 


(i) Using data from the VLA and VLB A we find the 
radio structure is from an AGN with a two-sided 
jet rather than being due to a dual AGN (Section 
3.1.1). 


(ii) Compared to the milliarcsecond absolute astrome¬ 
try provided by the VLBA the astrometry results 
suggest the hard X-ray, UV, optical, NIR, and ra¬ 
dio core emission are all coming from a single point 
source within < 0.2" (la. Section 3.1.3). We find 
no evidence of an offset between the radio and HST 



no evidence for a secondary AGN at 0.1" scales. 
This is consistent with ot her AO studies in the NIR 
( Imanishi fc Saito 2014) and mid-IR (8 — 13/am, 
Asmns et al.||2014p . 


(iii) With more than a factor of three deeper imaging in 
the hard X-rays with Chandra, we find the the pre¬ 
viously claimed dual AGN detection is most likely 
spurious resulting from the low number of X-ray 
counts (<160) at 6-7 keV and smoothing of data 
with a few counts per pixel on scales much smaller 
than the PSE (0.25" vs. 0.8" EWHM). We show 
that statistical noise in a single Chandra PSE gen¬ 
erates spurious dual peaks of the same separation 
(0.55 ± 0.07" vs. 0.6") and flux ratio (39 ± 9% vs. 
32% counts in the secondary) as the purported dual 
AGN (Section 3.1.4). 


(iv) We analyze 340 ks of recently obtained Chandra 
grating data in NGC 3393 finding that the pho- 
toionized emission dominates below 3 keV even in 
a 4" radius (Section 3.2.1). We detect a number of 
emission lines associated with strong photoionized 
emission at soft energies. We find large equivalent 
widths in the Ee Ka, Ee XXV a, Ee XXVI a con¬ 
sistent with a Compton-thick AGN. 


(v) We use NuSTAR spectroscopy to study the 
Compton-thick nature of NGC 3393, its torus ge¬ 
ometry, and its accretion rate (Section 3.2.2). We 
find both the BNtorus and MYtorus models sup¬ 
port a Compton-thick torus seen edge-on, with a 
large opening angle. Using the available NuSTAR, 
Chandra, XMM-Newton, Suzaku, BeppoSAXX-idy 
observations of NGC 3393, we find no evidence of 
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Figure 15. Tests for variability of N(^C 3393 in the X-rays. Upper Left. Relative X-ray flux normalized to the NuSTAR 2013 obser¬ 
vation, from various X-ray observatories from 2004-2013 for NGC 3393, calculated by fitting the data with the normalization for each 
telescope allowed to vary. The red points show normalized 3-8 keV flux while the blue points denote normalized 14-60 keV flux. We And the 
data show no variability below 10 keV between 2004-2011, followed by an increase of '^2 when the NuSTAR data was taken in 2013. Above 
10 keV, we And a similar increase of ^2 when comparing the stacked BAT 2004-2010 observations and NuSTAR. The lower sensitivity 
Suzaku PIN data shows no significant offset from NuSTAR. Bottom Left: Swift BAT 70-month 14-195 keV light curve spanning 2004-2010 
binned into 3 month intervals. We do not see any evidence of variability based on fitting with a constant value. Right: BNtorus model At 
with the different data sets. NuSTAR is shown in black (FPMA) and red (FPMB). Above 10 keV, Swift BAT (dark blue) and Suzaku PN 
(light green) are shown. Below 10 keV, XMM-Newton PN (cyan), MOS 1 (orange), and MOS 2 (pink), as well as Chandra data from 2004 
(green) and 2011 (pink) are shown. 


variability between 1997 and 2013 (Section 3.2.3). 
There is some weak evidence of variability above 
10 keV between the stacked Swift BAT observa¬ 
tions and NuSTAR at the l.Scr level dependent on 
the cross-calibration uncertainty. We find a high 
Eddington ratio of ?^0.2 based on the [Nev] and 
intrinsic 2-10 keV X-ray emission (Section 3.2.4). 

The most important factor in the likely spurious de¬ 
tection of a dual AGN in NGC 3393 was the low num¬ 
ber of X-ray counts (<160) at 6-7 keV and smoothing of 
data with a few counts per pixel on scales much smaller 
than the PSF (0.25" vs. 0.8" FWHM). Further Chandra 
searches for dual AGN with small angular separations 
(< 1") would benefit from studying sources with known 
optical or NIR counterparts such as found using high an¬ 
gular resolution HST or adaptive optics imaging. We are 
now extending this Chandra analysis to a larger sample 
of all Chandra observed Swift BAT detected AGN to test 
for any other spurious dual AGN detections, place upper 
limits on the frequency of closely separated dual AGN, 
and search for AGN that are significantly offset from the 
galaxy nucleus (Koss et ah, in prep). 

The combination of high-quality radio and X-ray data 
is useful for studying NGC 3393. The heavily inclined 
torus geometry inferred from the X-ray analysis agrees 


well with the fact that a disk water mega-maser has 
been detected in NGC 3393, since masing implies a 
nearl y edge-on view of the masing disk (Kondratko et al. 
2008). At a luminosity of 2 x lO^^ergs”^, we would 


expect a small opening angle of the torus due to the 
obsc ured fraction being n ear its peak at this luminosity 

(e-g. 

whet. 


Burlon et al. 2011). It is interesting to consider 
ler the large opening angle of the torus found in 
the X-rays is linked through feedback to the strong 
radio jets. Larger samples of radio bright obscured 
AGN will need to be studied to determine this link more 
thoroughly. 

NuSTAR has been critical for studies of local and dis- 


tant AGN above 10 keV (e.g., Balokovic et al. 2014 


Gandhi et al.||2014| Lansbnry et ai.||2014[ |Fiiccetti et al 


Arevalo et al.||2014| |Del Moro et al.||2014[ [Stern 


2014[ Arevalo et al.||2014| |Del Moro et al.||2014[ 

et al. 2014p. With the full coverage of the AGN from" 
0.2-100 keV with deep archival Chandra and NuSTAR 
observations we can assess the importance of this type 
of Compton-thick column density to different high red- 
shift X-ray AGN surveys. Previous studies have found 
that observations above 10 keV can be as successful at 
identifying Compton-thick AGN as deep X-ray obser¬ 
vations of these AGN probe softer X-ray energies (< 
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10 keV, e.g. Koss et ^^[2013 ). The intrinsic I/ 2 -iokeV = 
2.6 ± 0.3 X 10^^ erg s“^ is consistent with the median lu¬ 
minosity in the Chandra 4 Ms field at redshift z=1.6 
where much of black hole growth is thought to occur. As 
an unobscured source it could be detected up to redshift 
2^=5.0 at SN=3 if observed near the center of the field 
(< 3') where the PSF is sharpest and a 2" extraction 
radius encloses > 90% of the flux. We ran simulations of 
NGC 3393 observed on-axis with different satellites. In 
the Swift BAT all-sky survey, Compton-thick AGN like 
NGC 3393 are only detectable in the very nearby uni¬ 
verse {z <0.025), whereas with NuSTAR such objects 
are detectable out to z=0.075 with a 20 ks observation. 
This simulation suggests that near the center of the deep¬ 
est 4 Ms images, NGC 3393 would be detectable out to 
2:=0.55 at SN=3. However, the detection of the Fe K 
line at SN=3 and knowledge of the high level of obscu¬ 
ration and intrinsic luminosity would only be possible to 
redshifts 2 ; <0.2. Future observations with NuSTAR of a 
large sample of nearby Compton-thick AGN {z <0.5) are 
critical to understand how large a fraction of Compton- 
thick AGN might be hidden in the deepest Chandra im¬ 
ages. 
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